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Q Abstract 

C/3 The recent detection of both pressure and high-order gravity modes in the clas- 
, ^ , sical B-type pulsator 7 Pegasi offers promising prospects for probing its internal 
structure through seismic studies. To aid further modelling of this star, we 
^-H present the results of a detailed NLTE abundance analysis based on a large 
number of time-resolved, high-quality spectra. A chemical composition typical 
of nearby B-type stars is found. The hybrid nature of this star is consistent with 
ly-j its location in the overlapping region of the instability strips for j3 Cephei and 
slowly pulsating B stars computed using OP opacity tables, although OPAL cal- 
culations may also be compatible with the observations once the uncertainties 
in the stellar parameters and the current limitations of the stability calcula- 
T— I tions are taken into account. The two known frequencies /i = 6.58974 and 
, /2 = 0.68241 d^^ are detected in the spectroscopic time series. A mode iden- 
J>. tification is attempted for the low-frequency signal, which can be associated to 
a high-order y-mode. Finally, we re-assess the binary status of 7 Peg and find 
no evidence for variations that can be ascribed to orbital motion, contrary to 
previous claims in the literature. 



Individual Objects: 7 Pegasi. 



* Based on observations obtained at Vainu Bappu Observatory, Kavalur, 
India. 



2 



A spectroscopic study of the hybrid pulsator 7 Pegasi 



Introduction 

Of particular relevance for our understanding of the fundamental properties of 
stars on the upper main sequence are the so-called hybrid B pulsators, which 
simultaneously exhibit low-order pressure and high-order gravity modes charac- 
teristic of /? Cephei and slowly pulsating B stars (hereafter SPBs), respectively 
(in a similar vein, some A- and F-type main sequence stars also present 7 
Doradus- and S Scuti-like pulsations; e.g., Grigahcene et al. 2010). A num- 
ber of hybrid (3 Cephei/SPB pulsating candidates have been identified to date 
(e.g., De Cat et al. 2007, Pigulski Pojmariski 2008, Degroote et al. 2009), 
but their number is expected to grow dramatically in the future as intensive 
space observations with unprecedented duty cycle and photometric precision 
are being undertaken (e.g., Balona et al. 2011). The self-driven excitation of 
both pressure and gravity modes holds great asteroseismic potential because it 
offers an opportunity to probe both the stellar envelope and the deep internal 
layers. Such stars are hence prime targets for in-depth seismic modelling and, 
although challenging, their study promises to lead to significant advances in our 
understanding of the internal structure of main-sequence B stars (e.g., internal 
rotation profile, extent of the convective core; Thoul 2009). 

Of particular interest in this context is the bright 82 IV star 7 Peg (HR 
39, HD 886). Although for long considered as a classical /3 Cephei star, it has 
recently also been shown to exhibit gravity mode pulsations typical of SPBs, as 
first shown by Chapellier et al. (2006) using ground-based spectroscopic obser- 
vations. This result has recently been confirmed and extended by Handler et al. 
(2009) from high-precision photometry with the MOST satellite and a coordi- 
nated radial-velocity monitoring from the ground, with the detection of eight 
/3 Cep-like and six SPB-like oscillation modes (see also Handler 2009). With 
this peculiarity, 7 Peg enters the (so far) restricted club of hybrid pulsators and 
is expected to play in the future a pivotal role towards a better understanding 
of the physics of B stars. Modelling this star also has the potential to address 
even more fundamental issues, such as the reliability of opacity calculations 
(see the attempts in this direction by Walczak &i Daszynska-Daszkiewicz 2010 
or Zdravkov &i Pamyatnykh 2009). 

However, the reliability of the results provided by such seismic studies 
strongly hinges upon an accurate knowledge of both the fundamental parame- 
ters of the star under study (e.g., Teff) and its metal mixture. For this reason, 
an accurate determination of these quantities is of vital importance. In virtue of 
its brightness and initial status as a prototypical (3 Cephei star, many abundance 
works have been devoted to 7 Peg (e.g., the pioneering work of Aller 1949), 
but the analyses very often relied on non fully line-blanketed model atmospheres 
or fundamental parameters inferred from photometric calibrations and/or LTE 
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calculations (e.g., Gies &l Lambert 1992, Pintado & Adelman 1993, Ryans et 
al. 1996). These limiting assumptions cast some doubts on the reliability of 
these results. For instance, fitting of the wings of the Balmer lines using LTE 
synthetic spectra have been shown to systematically lead to an overestimation 
of the surface gravity (Nieva &l Przybilla 2007). Only a few studies have been 
conducted from an NLTE perspective (Andrievsky et al. 1999, Korotin et al. 
1999a, b, Morel et al. 2006), all but one of them only deriving the CNO abun- 
dances.^ Here we present a detailed NLTE abundance analysis of this star in 
an effort to re-examine its fundamental properties and its position relative to 
the theoretical instability domains for SPBs and /3 Cephei stars. We also take 
advantage of the large number of time-resolved, high-quality spectra collected 
to re-assess its binarity and to attempt to identify the most prominent modes 
visible in spectroscopy. 



Observations and data reduction 

The spectroscopic observations were obtained at Vainu Bappu Observatory 
(VBO) located in Kavalur (India) using the fiber-fed echelle spectrograph at- 
tached to the prime focus of the 2.3-m Vainu Bappu telescope (VBT). The 
wavelength range was 4000-8000 A for the 2 K x 4 K CCD detector (spread 
over 45 orders) and 4200-7000 A for the 1 K x 1 K CCD detector (spread over 
nearly 25 orders with the echelle gaps). The resolving power estimated from the 
arc spectra is i? = A/AA w 60 000. Further details regarding the instrumental 
set up are given in Rao et al. (2005). 

A total of 163 spectra were collected from September 2007 to January 
2009 during various observing runs. The journal of observations is presented 
in Table 1. The exposure time ranged from 10 to 15 minutes depending on 
the sky conditions (i.e., clear sky or sky with thin clouds) and the signal-to- 
noise ratio was generally above 200. The data were reduced and analysed 
using standard IRAF (Image Reduction and Analysis Facility)^ tasks. The basic 
steps of the data reduction included trimming, bias frame subtraction, scattered 
light removal, flat fielding, order extraction, wavelength calibration (using ThAr 
lamps) and finally continuum rectification. 



'^Note that the NLTE abundances of Gies & Lambert (1992) are not based on full 
NLTE linc-forination calculations, but arc derived instead from the LTE values assuming 
theoretical NLTE corrections appropriate for a star with parameters typical of 7 Peg. 

"^IRAF is distributed by the National Optical Astronomy Observatories, operated by the 
Association of Universities for Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation. 
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Table 1: Journal of observations. 



Civil Date 


Number of spectra 


Civil Date 


Number of spectra 


20 Sep. 


07 


10 


02 Feb. 08 


4 


28 Sep. 


07 


19 


03 Jun. 08 


4 


02 Oct. 


07 


24 


04 Jun. 08 


4 


04 Oct. 


07 


11 


05 Jun. 08 


2 


05 Oct. 


07 


28 


24 Aug. 08 


5 


01 Jan. 


08 


3 


31 Aug. 08 


3 


02 Jan. 


08 


11 


04 Oct. 08 


7 


03 Jan. 


08 


5 


05 Oct. 08 


11 


16 Jan. 


08 


3 


13 Jan. 09 


3 


17 Jan. 


08 


3 


14 Jan. 09 


3 




Total 


163 



Line-profile variations 

It should be noted that because of different instrumental settings, not all spec- 
tral lines were systematically covered during the observations. We therefore 
focus in the following on the strong spectral lines that were the most exten- 
sively observed, namely the Si III triplet between 4552 and 4575 A and the C 
II AA5143, 5145 doublet, with 80 and 131 exposures, respectively. 

Binarity 

Until recently, 7 Peg was believed to be a spectroscopic binary whose orbital 
period was, however, disputed (Porb = 370.5 days, Chapellier et al. 2006; 
Porb = 6.816 days, Harmanec et al. 1979, Butkovskaya & Plachinda 2007). 
Contrary to previous authors. Handler et al. (2009) suggested 7 Peg to be 
a single star and explained the claimed orbital variations as due to 5-mode 
pulsation. 

To re-assess the possible binarity of this object, we compared our radial 
velocities (RVs) of the C II A5145 line with the orbital solution proposed by 
Chapellier et al. (2006). Our heliocentrically-corrected RVs vary between —1.5 
and 7 kms~^ (with a typical accuracy on the individual measurement of ~0.3 
kms~^) mainly as a result of the well-known, dominant radial mode with a 
period of about 0.15 day (see below). There is no evidence for any long- 
term trends, despite the fact that our observations cover some phase intervals 
(most notably around 0~O.93) where rapid variations amounting to up to ~12 
kms~^ are expected according to the ephemeris of Chapellier et al. (2006). 
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This rules out the possibility that 7 Peg is an eccentric binary with an orbital 
period of about one year, confirming the conclusion of Handler et al. (2009). 
On the other hand, these authors favoured the pulsation interpretation for the 
~6.8 days period since they noticed that the frequency 1/6.816 d^^ with an 
amplitude of ~0.8 kms~^ (Butkovskaya &i Plachinda 2007), is a one-day alias 
of the frequency 0.8533 d~^, which lies well within the domain of 5-mode 
frequencies. Finally, our data also do not support the large, abrupt RV changes 
reported by Butkovskaya it Plachinda (2007). 

Frequency analysis and mode identification 

To further investigate the line-profile variability in 7 Peg, we used the software 
package FAMIAS^ (Zima 2008). We computed the first three velocity moments 
< >, <v'^ > and <v'^ > (see Aerts et al. 1992 for a definition of the 
moments of a line profile) of the Si III triplet and the C II doublet with the 
aim of performing a frequency analysis. The integration limits for computing 
the moments were dynamically chosen by sigma clipping to avoid the noisy 
continuum. 

In the first moment < >, which is the RV placed at average zero, we 
detected the well-known dominant radial mode, which can readily be seen in 
the line-profile variations affecting the Si III lines (Fig. 1), followed by a one- 
day alias of the known highest amplitude low frequency mode. Both pulsation 
modes were first discovered by Chapellier et al. (2006) and afterwards con- 
firmed by Handler et al. (2009), who definitely proved the hybrid (3 Cep/SPB 
nature of 7 Peg. Since the frequency values of Handler et al. (2009) are more 
accurate than ours, we adopted them in our study, i.e., /i = 6.58974 and 
/2 = 0.68241 d~^. The corresponding amplitudes in our dataset are 3.39 and 
0.69 kms^^ with an error of 0.07 kms^^, for /i and /2, respectively. No addi- 
tional frequencies could be detected according to the 4 signal-to-noise criterion 
of Breger et al. (1993). The latter was tested by computing the noise level 
for different box intervals (between 1 and 10 d~^) centred on the considered 
frequency. 

In an attempt to detect other periodicities, we performed a frequency search 
on the spectra by means of a two-dimensional Fourier analysis. In this way, we 
again detected /i and /2 but no additional frequencies. A non-linear multi- 
periodic leastsquares fit of a sum of sinusoidals is then computed with the 

Levenberg-Marquardt algorithm (adapted from Press et al. 2007). This fitting 
is applied for every bin of the spectrum separately according to the formula 
Z + ^ . Ai sin[27r(/it -I- </>,)] , where Z is the zeropoint, and Ai, fi, and 0i are 



^FAMIAS has been developed in the framework of the FP6 European Coordination 
Action HELAS - http://www.helas-eu.org/ 
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Figure 1: Examples of line-profile variations affecting Si III A4567.8, 4574.7 and 4813.3 
on 5 October 2007. 

the amplitude, frequency, and phase of the z-th frequency, respectively. The 
amplitude and phase distributions across the C II A5145 A line are shown in 
Fig. 2. 

To identify the modes associated to /i and /2, we used the Fourier pa- 
rameter fit method (FPF method; Zima 2006). The wavenumbers {£,m) and 
other continuous parameters are determined in such a way that the theoretical 
zeropoint, amplitude and phase values across the profile best fit the observed 
values. The fitting is carried out by applying genetic optimization routines in 
a large parameter space. From multicolour photometric time series. Handler 
(2009) unambiguously identified /i as a radial mode while two possibilities 
remain for /2, which is a £2 = 1 or a £2 = 2 mode. Our spectroscopy corrob- 
orates that (£i,mi) — (0,0), as illustrated in Fig. 2. To identify the values of 
(£2,^2), we only allowed £2 to be 1 or 2. Such an approach, which consists 
in adopting the £-\/a\ues obtained from photometry, already proved to be suc- 
cessful for spectroscopic mode identification of other /? Cephei stars (Briquet 
et al. 2005; Mazumdar et al. 2006; Desmet et al. 2009). In Table 2 the best 
parameter combinations are listed for each couple {£2, ^2). The smaller the 
value, the better the solution. The most probable and equally good solutions 
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V(kms'') V(kms"^) 



Figure 2: Amplitude and phase distributions (points with error bars) across the C II 
A5145 A line, for the frequencies /i = 6.58974 and /2 = 0.68241 d"\ The bestfit 
model is represented by full lines. For /i, the mode is radial. For /2, the differences 
between the best solutions listed in Table 2 for (£2, ?n2)=(l,l), (2,-1) and (2,-2) 
are smaller than the thickness of the lines. The amplitudes are expressed in units of 
continuum and the phases in tt radians. 

are (^'2,^2) — (1.1). (2,-2) or (2,-1) (where a positive m-value denotes a 
prograde mode). The best models for the amplitude and phase across the line 
profile are illustrated in Fig. 2. A by-product of the mode identification is the 
derivation of the equatorial rotational velocity Woq. By using the x^-values as 
weights, we constructed a histogram for (see Fig. 3), as in Desmet et al. 
(2009). We computed it by considering the solutions with (^2,^2) = (1,1), 
(2,-2) and (2,-1). By calculating a weighted mean and standard deviation, 
we obtained Weq = 6±1 kms~^ (a much less probable value is ~12 kms~^). 

Abundance analysis 
Analysis tools 

The atmospheric model calculations were performed under the assumption of 
LTE, whereas a full NLTE treatment was adopted for the line formation. Such 
a hybrid approach has been shown to be adequate for early B-type stars on the 
main sequence (Nieva h. Przybilla 2007). First, the ATLAS9 code^ (Kurucz 
1993) is used to compute the hydrostatic, plane-parallel and fully line-blanketed 



*http: / /kurucz. harvajd.edu/ 
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Table 2: Mode parameters derived from the FPF mono-mode method for /2. as is 
the surface velocity amplitude (km s~^); i is the stellar inclination angle in degrees; 
vsini is the projected rotational velocity, a is the width of the intrinsic profile, both 
expressed in km s~^. 





(^2,^2) 


as 


i 


vsini 


a 


3.59 


(1,1) 


3.9 


14 


1.3 


6.9 


3.60 


(2,-2) 


2.5 


25 


2.6 


6.5 


3.62 


(2,-1) 


2.5 


84 


6.2 


6.1 


4.77 


(1,0) 


6.2 


86 


5.8 


5.3 


4.89 


(0,0) 


2.8 




1.0 


6.9 


4.91 


(2,0) 


0.4 


26 


4.7 


6.6 


5.29 


(2,1) 


5.1 


87 


1.0 


6.4 


5.30 


(1,-1) 


0.9 


63 


1.0 


7.1 


5.84 


(2,2) 


7.4 


15 


1.0 


6.1 



35 




Vg (km s"'') 



Figure 3: Histogram for the equatorial rotational velocity derived from the FPF 
method. 
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LTE atmospheric models. Grids with the new opacity distribution functions 
(ODFs) assuming the solar abundances of Grevesse &i Sauval (1998) have been 
used (Castelli & Kurucz 2004). 

Further, in order to obtain the NLTE abundances we made use of the lat- 
est version of the NLTE line-formation codes DETAIL/SURFACE (Butler & 
Giddings 1985, Giddings 1981). DETAIL provides the solution of the radiative 
transfer and statistical equilibrium equations, while the emergent spectrum is 
calculated by SURFACE. The line atomic data are taken from the NIST and 
VALD databases. Care has been taken to only retain features that are unblended 
In the relevant temperature range. 

Determination of the atmospheric parameters 

Our abundance analysis is based on the average of a large number of time- 
resolved spectra (selected to have the highest S/N) to ensure that the param- 
eters and abundances we derive are representative of the values averaged over 
the whole pulsation cycle. 

The effective temperature, Tgg, was estimated from the Si ll/lll ionisation 
balance. We cannot rely upon Si lll/IV ionisation balance, as no Si IV lines 
are visible. On the other hand, in this Toff range no other chemical elements 
have lines of two adjacent ionisation stages that can be measured. Numerous 
Fe III lines and two very weak Fe II lines (Fe II A5018 and Fe II A5169) with an 
EW in the range 10-15 mA are present In our spectra. However, constraining 
Tefj from ionisation balance of iron was not possible owing to the rudimentary 
nature of our Fe II model atom, which only includes 8 levels. In contrast, the 
detailed treatment of the Fe III ion (264 levels) allows us to reliably estimate 
the iron abundance from the analysis of the Fe III features. Figure 4 shows 
examples of calibrations between the Si line ratios and Teff. Ten line ratios 
have been used for the temperature estimation and the dispersion between the 
various Tefi values obtained was considered as representative of the uncertainty 
in this parameter. 

Unfortunately, we were unable to estimate the surface gravity, log 51, from 
the fitting of the collisionally-broadened wings of the Balmer lines, as no profiles 
were completely covered because of echelle gaps. Therefore, this parameter 
was estimated from Stromgren uvby/3 photometry based on the calibrations of 
Castelli (1991). The ci and (3 colour indices were taken from the Catalogue 
of uvbyP data of Hauck 8i Mermilliod (1998). The observed ci colour index 
was dereddened to obtain cq by means of empirical calibrations (details can 
be found in Castelli 1991 and references therein). Slightly different sets of 
colour indices (ci, /3) are available in the literature (Crawford et al. 1971, 
Sterken & Jerzykiewicz 1993). The surface gravity was computed using all the 



10 A spectroscopic study of the hybrid pulsator 7 Pegasi 



9:8 



5056/4568 

t = Q5 

1=10 

, ■"■"■|=20 


\5056/4574 

\ 


5056/4813 


^56/4829 

"\ 

\ 
\ 

,v 


\5056/5739 
\ 


6371 /4568 

V 


6371 /4574 

\ 


6^71/4813 

\ 


^1^71/4829 

\ 

\ 


1^ 6371/5739 



20 22 24 26 26 22 24 26 26 22 24 26 26 22 24 26 26 22 24 26 28 
Teff{kK) Teff(kK) Teff(kK) Teff(kK) Teff(kK) 



Figure 4: Examples of calibrations between various Si line ratios and the effective 
temperature, as a function of the microturbulence ^ (red lines: ^ = 5 km s^^; light 
green lines: ^ = 10 km s^^; blue lines: ^ = 15 km s^^; dark green lines: ^ = 20 km 
s^^). We have adopted logg = 3.7 and ^ = 5 km s^^, as appropriate for this star 
(see text). 



possible combinations and the scatter was taken as the typical error bar. The 
simultaneous estimation of the effective temperature and surface gravity, along 
with the error limits, is illustrated in Fig. 5. 

Another important parameter which is required for the abundance analysis 
is the microturbulent velocity, ^, which was determined from the II features 
requiring that the abundances are independent of the line strength (Fig. 6). The 
uncertainty in the microturbulence was estimated by varying this parameter until 
the slope of the log e(0)-log e(EW/A) differs from zero at the 3ct level. 

Finally, the total amount of line broadening, vt, is determined from line- 
profile fitting of a number of isolated metallic lines (instrumental broadening was 
estimated from the calibration lamps). Although rotational broadening clearly 
dominates over pulsational broadening in this star whose dominant pulsation 
mode is radial, the fact that this estimate {vt = 8±1 km s^^) is marginally 
higher than the equatorial rotational velocity inferred from the moment analysis 
(ueq = 6±1 km s^^) suggests that broadening arising from non-radial pulsations 
may not be completely negligible. A comparison between the observed and fitted 
profiles is presented in the Appendix. 

The atmospheric parameters, along with their uncertainties, are presented in 
Table 3. A detailed comparison with previous results in the literature is deferred 
to the last section. 
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Figure 5: Tcff-logij plane for the simultaneous determination of the effective tem- 
perature and surface gravity. The group of nearly vertical lines represents the loci 
satisfying the Si I I/Si III ionisation equilibria, while the thick solid line represents the 
log(? determined from Stromgren photometry. The middle point of the intersection 
of the two sets of lines simultaneously provides the value of Toft and log g. 



Determination of the elemental abundances 

The NLTE abundances of He, C, N, 0, Ne, Mg, Al, Si and Fe have been 
determined from a classical curve-of-growth analysis using the EWs measured 
through direct integration. The mean abundances are presented in Table 4, 
while the results for each transition considered are given in the Appendix (along 
with the EW measurements). To estimate the abundance uncertainties, we have 
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Figure 6: Estimation of the microturbulent velocity from the II lines. 



first calculated the errors arising from the uncertainties in the atmospheric pa- 
rameters (i.e., (TTcff , ciogg and cr^). They were derived by computing the abun- 
dances using models with atmospheric parameters deviating from the nominal 
values by the uncertainties tabulated in Table 3. We also considered the fact 
that the determinations of Tcs and logy are strongly coupled. Finally, the total 
uncertainty, Ctotah is obtained by quadratically summing up these errors and 
the line-to-line scatter ((Tint)- 

After completion of the analysis, we became aware of a study (Simon-Dfaz 
2010) reporting problems with the modelling of some silicon lines using the 
model atom implemented in the NLTE code FASTWIND (Puis et al. 2005). 
Because our model atom is similar in many respects (Morel et al. 2006), 
we have examined to what extent this would affect our determination of the 
atmospheric parameters and ultimately abundances by repeating the abundance 
analysis excluding the spectral lines that may not be properly modelled (namely 
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Table 3: Derived atmospheric parameters for 7 Peg. 



Teff (K) 
logTeff 

log .9 

^ (km s-i) 
Vt (km s~^) 



22650±650 
4.355±0.013 
3.73±0.08 

8±1 



Table 4: Mean NLTE abundances (the number of lines used for each element is given 
in brackets) with the usual notation assuming loge(H)=12, and details on the error 
budget. The last column gives the total abundance uncertainty. 



Mean Error estimation 



abundance 


(Tint 


0"Tcff 


criogg 




""Tcff/logg 


Ctotal 


loge(He) 


10.92 


(8) 


0.24 


0.04 


0.02 


0.05 


0.03 


0.25 


loge(C) 


8.26 


(9) 


0.05 


0.00 


0.00 


0.01 


0.01 


0.05 


loge(N) 


7.62 


(22) 


0.08 


0.07 


0.03 


0.03 


0.03 


0.12 


loge(O) 


8.43 


(23) 


0.05 


0.15 


0.06 


0.06 


0.09 


0.20 


loge(Ne) 


8.26 


(8) 


0.05 


0.05 


0.01 


0.01 


0.04 


0.08 


loge(Mg) 


7.60 


(1) 




0.06 


0.01 


0.16 


0.06 


0.18 


loge(Al) 


6.30 


(3) 


0.00 


0.04 


0.00 


0.09 


0.02 


0.10 


loge(Si) 


7.14 


(8) 


0.18 


0.06 


0.02 


0.09 


0.03 


0.21 


loge(Fe) 


7.30 


(24) 


0.09 


0.08 


0.05 


0.07 


0.03 


0.15 



Si II A5056 and Si III A4813, 4819, 4829). This leads to a slightly lower 
temperature (see Fig. 5) and surface gravity: Tdi — 22300 K and log g — 
3.72 (the microturbulence remains unchanged). The abundances determined 
using these parameters are provided in Table 5. The differences with the mean 
abundances derived assuming the default parameters (Table 4) remain within 
0.1 dex, except in the case of Si where it amounts to about 0.18 dex. 



Discussion 
Chemical composition 

Figure 7 shows a comparison between our abundances (Table 4) and previous 
NLTE results in the literature (Gies & Lambert 1992, Andrievsky et al. 1999, 
Korotin et al. 1999a, b. Morel et al. 2006, Morel & Butler 2008), as well 
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Table 5: Mean NLTE abundances (the number of lines used for each element is given 
in brackets) assuming T^g = 22300 K, logs' = 3.72 and ^ = 1 km s"^ 



Mean 

abundance 



loge(Hc) 


10.88 


(8) 


loge(C) 


8.25 


(9) 


loge(N) 


7.65 


(22) 


loge(O) 


8.52 


(23) 


loge(Ne) 


8.24 


(8) 


loge(Mg) 


7.56 


(1) 


loge(Al) 


6.32 


(3) 


loge(Si) 


7.32 


(4) 


log e(Fe) 


7.30 


(24) 



as the new solar abundances of Asplund et al. (2009). There is an overall 
satisfactory agreement with these previous studies. However, all these results 
taken at face value would suggest that 7 Peg is metal poor with respect to the 
Sun, a fact which is obviously contrary to the expectations for a young star 
in the solar neighbourhood. This is likely connected to the long-standing and 
more general problem affecting most abundance analyses of OB stars, which 
yield abundances significantly lower than solar (e.g., Morel 2009 for a review). 
Although the origin of this problem remains unclear, improvements in the atomic 
data and determination of the atmospheric parameters have been claimed to 
solve most of the discrepancy (Przybilla et al. 2008, Simon-Diaz 2010). In 
view of this rather unsatisfactory situation, a sound assumption may be to use 
in further theoretical modelling the abundances determined for a small sample 
of nearby B-type stars by Przybilla et al. (2008) or the solar abundances of 
Asplund et al. (2009). 

However, a robust result is that the chemical composition of 7 Peg does 
not significantly differ from that of early B-type stars analysed using similar 
techniques, whether they are known as pulsating or not (Morel et al. 2008). 
The existence of microscopic diffusion processes in hybrid pulsators leading to 
an accumulation of iron in the driving zone has often been invoked to explain 
the unexpected excitation of some specific pulsation modes (e.g., Pamyatnykh 
et al. 2004), but the lack of any abundance peculiarities does not seem at 
first glance to support this possibility unless one supposes that the changes 
at the surface remain below the detection limits. On the other hand, the 
[N/0] and [N/C] logarithmic abundance ratios (-0.81±0.24 and -0.64±0.13 
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Figure 7: Comparison between our abundances and previous NLTE results in the 
literature (Gies & Lambert 1992, Morel et al. 2006, Morel &i Butler 2008, Andrievsky 
et al. 1999, Korotin et al. 1999a, b). The data points dubbed 'Odessa group' refer 
to the results of Andrievsky et al. (1999) and Korotin et al. (1999a, b). The solar 
abundances of Asplund et al. (2009) are indicated with different symbols. 



dex, respectively) are identical to within the errors to the solar values (Asplund 
et al. 2009). Our study therefore confirms that there is no evidence in 7 Peg 
for CNO-processed material dredged up to the surface because of deep mixing. 
The significant boron depletion (Proffitt &l Quigley 2001) suggests, however, 
that shallow mixing is already taking place in the superficial layers. 



Comparison with theoretical instability strips 

As 7 Peg is one of the rare stars known to date to present both high-order g 
modes and low-order p and g modes, it is of interest to examine whether this 
property is consistent with the theoretical expectations. The position of 7 Peg 
in the (logTcff , log g) plane is shown in Fig. 8, along with the instability domains 
for (3 Cephei and SPB-like pulsation modes computed with the solar abundances 
of Asplund et al. (2005) and for two different opacity tables, OP and OPAL 
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(see Miglio et al. 2007 for details). At face value, this figure suggests that 
the hybrid nature of the pulsations in 7 Peg is more consistent with instability 
strips computed with OP opacities, as the OPAL ones do not predict this star 
to exibit SPB-like pulsations. 



GO 

tsfl 
O 



3.5 



4.5 




4.4 



4.2 
log(Teff) 



Figure 8: Position of 7 Peg in the (logTcff, log 3) plane (green filled circle), along 
with the instability strips for p Cephei- (red lines) and SPB-like (blue lines) pulsation 
modes {I < 4) computed using the solar abundances of Asplund et al. (2005) and 
the OP and OPAL opacity tables (Miglio et al. 2007). The neon abundance is taken 
from Cunha et al. (2006). The open circles indicate results for 7 Peg in the literature 
(Fitzpatrick & Massa 2005, Gies &i Lambert 1992, Korotin et al. 1999a, Martin 2004, 
Morel et al. 2006 [based on Si ionisation equilibrium]. Morel & Butler 2008 [based 
on Ne ionisation equilibrium], Niemczura &i Potubek 2006, Pintado &i Adelman 1993, 
Ryans et al. 1996). The two open squares refer to the results for two other hybrid 
pulsators, v Eri and 12 Lac (Morel et al. 2006). 



As can also be seen in this figure, the same conclusion holds for the two 
other best-studied hybrid pulsators, v Eri (Jerzykiewicz et al. 2005) and 12 Lac 
(Handler et al. 2006), that have been analysed using similar techniques as in 
the present work (Morel et al. 2006). Identical parameters within the errors 
were obtained for 1/ Eri by De Ridder et al. (2004) from a full spectroscopic 
analysis. Lefever et al. (2010) analysed 7 Peg, ly Eri and 12 Lac using the 
same data as Morel et al. (2006), but with a different code, and also found 
very similar results, except a temperature lower by 1500 K in 12 Lac. These 
two studies also determined T^s and log g from Si ionisation balance and fitting 
of the Balmer lines. On the other hand, and although they may be regarded as 
less reliable as they heavily rely on only one very weak Ne II line, slightly lower 
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Teff values were obtained for these three stars by Morel &i Butler (2008) based 
on Ne ionisation equilibrium. Finally, Gies &i Lambert (1992) found significantly 
higher Teff values for u Eri and 12 Lac based on photometric data, but it is 
likely that their temperature scale is too hot (see discussion in Lyubimkov et al. 
2002). This also leads to surface gravities that are higher than other results in 
the literature. 

While the opacity values prove critical for the excitation, or lack thereof, of 
oscillations in B stars, the location of the SPB instability domain is relatively 
insensitive to the chemical mixture used (Miglio et al. 2007). Using the recently 
revised abundances of Asplund et al. (2009) is therefore not expected to change 
this picture significantly (an identical situation is indeed encountered when using 
the more different solar mixture of Grevesse &l Noels 1993). 

To investigate to what extent our conclusions regarding the ability of the 
opacity calculations to reproduce the observations is robust against the choice 
of logTeff and logg, we show in Fig. 8 other results for 7 Peg from the liter- 
ature. Although a variety of techniques have been used for the determination 
of the parameters, almost all these results rely on photometric data and/or 
LTE methods. It should also be noted that these estimates are not completely 
independent, as the same photometric calibrations (and data) have often been 
employed. Despite these caveats, the gravity seems well constrained to around 
log 5 ~ 3.8 if one excludes the discrepant value of Pintado ii Adelman (1993) 
based on LTE fitting of the H7 line. The effective temperatures are gener- 
ally based on narrow-band photometric data and appear slightly lower than our 
value based on Si ionisation equilibrium (keeping in mind that, as discussed 
above, rejecting some potentially problematic Si lines would lower our T^g by 
about 350 K). Taking into account the still relatively large inaccuracy of the 
Teff value, we conclude, solely based on the position of 7 Peg with respect to 
the instability strips, that the OP opacities are only marginally preferred to the 
OPAL ones. Accurate positions for a sufficiently large sample of stars may lead 
to stronger statements. It should also be kept in mind, however, that some 
important physical processes that are not currently implemented in the sta- 
bility analysis may substantially change the results of the models (e.g., stellar 
rotation; Townsend 2005). 

To conclude, although it is clear that several theoretical aspects related to 
the excitation/damping of oscillations in B stars may greatly benefit from the 
constraints offered by the further identification of hybrid pulsators and from a 
detailed study of their pulsation properties, an accurate determination of their 
fundamental parameters remains essential for further progress. 
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1. Appendix 



He I He I He I He I He I He I He I He I 




4387.929 4437.551 4471.480 4713.15 5015.678 5047.738 5875.62 6678.154 



Figure 9a: Comparison between the observed and synthetic profiles (represented by 
the dots and full lines, respectively) computed for the abundance yielded by the cor- 
responding He I line. The synthetic spectra have been convolved with a rotational 
broadening function with vt=8 km s^^. 



C II C II C II C II C II C II C II C II C II 




5133.11 5137.26 5139.17 5143.40 5145.16 6783.91 6787.21 6791.47 6800.69 



Figure 9b: Same as Fig. 9a, but for the C II lines. 
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Figure 9c: Same as Fig. 9a, but for the N II lines. 
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II II II II II II II II 




4185.44 4319.78 4366.70 4414.90 4416.97 4452.38 4590.97 4596.07 
II II II II II II II II 




4641.81 4649.13 4661.63 4673.73 4676.23 4696.32 4701.44 4705.35 
II II II II II II II 




4890.85 4906.83 4941.07 4943.00 4955.71 5190.49 5206.64 



Figure 9d: Same as Fig. 9a, but for the II lines. 



Ne I Ne I Ne I Ne I Ne I Ne I Ne I Ne I 




0.8 [ \ I \ I \ I \ I \ I \ I \ I \ 

5852.488 6143.063 6163.594 6217.287 6334.428 6382.991 6402.246 6506.528 



Figure 9e: Same as Fig. 9a, but for the Ne I lines. 
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Mg II Al III Al III Al III 




44S1.20 4479, S6 4512-54 4529, 04 



Figure 9f: Same as Fig. 9a, but for the Mg II (left) and Al III (right) lines. 



Si II Si II Si III Si III Si III Si III Si III Si III 



5056.15 5371.371 4567.84 4574.757 4813.33 4819.77 4829.073 5739.734 



Figure 9g: Same as Fig. 9a, but for the Si II and Si III lines. 
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Figure 9h: Same as Fig. 9a, but for the Fe III lines. 
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Table 6: EW measurements and abundance results for each transition, along with the 
error estimates. 



Ion 


Transition 

(A) 


EW 
(mA) 


log e 
(dex) 




Error estimation 

O'logg CTf CTxeff/logg 




He I 


4387.93 


820 


11.29 


0.08 


U.Uo 


U.UU 


U.U4 


0.09 




4437.55 


109 


10.98 


0.06 


n on 
U.UU 


U.UD 


U.Uo 


0.10 




4471.48 


1255 


11.09 


0.05 


0.04 


0.00 


0.00 


0.06 




4713.15 


270 


10.92 


0.07 


0.00 


0.07 


0.07 


0.12 




5015.68 


274 


10.76 


0.00 


n no 

u.uy 


n 1 n 
U.iU 


n nn 
U.UU 


0.13 




5047.74 


170 


11.04 


0.06 


U.UU 


n nfi 
U.Uo 


n n*^ 

U.UD 


0.10 




5875.62 


605 


10.53 


0.00 


U.UU 


n nn 
U.UU 


n nn 
U.UU 


0.00 




6678.15 


560 


10.76 


0.00 


U.UU 


n 1 n 

U.iU 


n nn 

U.UU 


0.10 


C II 


5133.11 


62 


8.22 


0.01 


U.UU 


n ni 
U.Ui 


n ni 
U.Ui 


0.02 




5137.26 


9 


8.31 


0.00 


n nn 

U.UU 


n nn 

U.UU 


n ni 

U.Ui 


0.01 




5139.17 


14 


8.33 


0.01 


n ni 
U.Ui 


n m 
U.Ui 


n no 
U.Uz 


0.03 




5143.40 


33 


8.26 


0.01 


n nn 
U.UU 


n ni 
U.Ui 


n ni 
U.Ui 


0.02 




5145.16 


58 


8.25 


0.01 


n ni 
U.Ui 


n OQ 
U.Uo 


n ni 
U.Ui 


0.03 




6783.91 


50 


8.25 


0.00 


n ni 
U.Ui 


n nQ 
U.Uo 


n ni 
U.Ui 


0.03 




6787.21 


17 


8.30 


0.01 


n nn 

U.UU 


n ni 

U.Ui 


n ni 

U.Ui 


0.02 




6791.47 


18 


8.23 


0.00 


n ni 

U.Ui 


n ni 

U.Ui 


n ni 

U.Ui 


0.02 




6800.69 


15 


8.20 


0.01 


n nn 
U.UU 


n m 
U.Ui 


n ni 
U.Ui 


0.02 


N II 


4227.74 


16 


7.53 


0.05 


n ni 
U.Ui 


n no 
U.Uz 


n nQ 
U.Uo 


0.06 




4236.99 


38 


7.48 


0.06 


n no 
U.Uz 


n m 
U.Ui 


n nQ 
U.Uo 


0.07 




4607.16 


40 


7.73 


0.07 


n no 

U.U/ 


n nQ 
U.Uo 


n nQ 
U.Uo 


0.08 




4643.09 


42 


7.62 


0.06 


n nQ 
U.Uo 


n r\A 
U.U4 


n nQ 
U.Uo 


0.08 




4779.72 


11 


7.66 


0.05 


n ni 
U.Ui 


n nn 
U.UU 


n no 
U.Uz 


0.05 




4788.14 


18 


7.70 


0.06 


n no 
U.Uz 


n no 
U.Uz 


n nQ 
U.Uo 


0.07 




4987.38 


10 


7.70 


0.05 


n 0? 


ni 


n 0*^ 


0.06 




4994.36 


19 


7.53 


0.07 


0.02 


0.02 


0.04 


0.09 




5001.30 


83 


7.59 


0.09 


0.03 


0.05 


0.05 


0.12 




5002.70 


16 


7.72 


0.05 


0.03 


0.02 


0.02 


0.06 




5005.17 


53 


7.55 


0.10 


0.03 


0.06 


0.06 


0.13 




5007.33 


28 


7.56 


0.08 


0.02 


0.03 


0.04 


0.10 




5010.62 


29 


7.65 


0.06 


0.03 


0.03 


0.02 


0.08 




5025.66 


13 


7.76 


0.07 


0.02 


0.01 


0.04 


0.08 




5045.10 


39 


7.65 


0.06 


0.04 


0.04 


0.03 


0.09 




5495.65 


10 


7.53 


0.07 


0.02 


0.01 


0.04 


0.08 




5666.63 


50 


7.57 


0.09 


0.04 


0.05 


0.05 


0.12 




5676.01 


35 


7.64 


0.08 


0.04 


0.04 


0.04 


0.11 




5679.55 


71 


7.58 


0.10 


0.04 


0.07 


0.05 


0.14 




5686.21 


21 


7.50 


0.07 


0.03 


0.02 


0.03 


0.08 




5710.77 


33 


7.76 


0.08 


0.04 


0.04 


0.04 


0.11 




5747.30 


8 


7.59 


0.06 


0.03 


0.01 


0.03 


0.07 
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II 


4185.44 


30 


8.42 





15 


0.05 


0.07 


0.10 


0.20 




4319.78 


60 


8.60 





15 


0.06 


0.08 


0.10 


0.21 




4366.70 


48 


8.40 





15 


0.06 


0.08 


0.09 


0.20 




4414.90 


70 


8.37 





15 


0.05 


0.09 


0.10 


0.21 




4416.97 


56 


8.41 





15 


0.06 


0.09 


0.09 


0.21 




4452.38 


25 


8.44 





13 


0.05 


0.04 


0.08 


0.17 




4590.97 


52 


8.48 





18 


0.06 


0.08 


0.12 


0.24 




4596.07 


50 


8.47 





16 


0.06 


0.05 


0.10 


0.20 




4641.81 


68 


8.40 





18 


0.06 


0.10 


0.11 


0.24 




4649.13 


92 


8.51 





17 


0.07 


0.14 


0.11 


0.26 




4661.63 


47 


8.37 





16 


0.06 


0.08 


0.10 


0.21 




4673.73 


17 


8.42 





12 


0.05 


0.03 


0.08 


0.16 




40 / D.ZO 


A K 
40 


Q Afi 


u 


1 

10 


U.UO 


n n7 
u.u / 


U.iU 


n 91 




4696.32 


10 


8.38 





12 


0.05 


0.02 


0.07 


0.15 




4701.44 


10 


8.39 





14 


0.04 


0.02 


0.09 


0.17 




4705.35 


40 


8.38 





16 


0.05 


0.07 


0.11 


0.21 




4890.85 


12 


8.44 





14 


0.06 


0.02 


0.08 


0.17 




4906.83 


19 


8.42 





16 


0.06 


0.04 


0.09 


0.20 




4941.07 


19 


8.43 





16 


0.05 


0.04 


0.10 


0.20 




4943.00 


27 


8.41 





17 


0.05 


0.06 


0.11 


0.22 




4955.71 


8 


8.42 





13 


0.05 


0.02 


0.08 


0.16 




CI on AO 





Si AO 

o.4y 





14 


U.UO 


u.u± 


nsi 

U.UO 


17 
U. ± / 




5206.64 


16 


8.37 





15 


0.05 


0.03 


0.10 
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